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Abstract. Scalable quantum networks require the capability to create, store and distribute entanglement 
among distant nodes (atoms, trapped ions, charge and spin qubits built on quantum dots, etc.) by means of 
photonic channels. We show how the entanglement between qubits and electromagnetic field modes allows 
generation of entangled states of remotely located qubits. We present analytical calculations of linear 
entropy and the density matrix for the entangled qubits for the system described by the Jaynes-Cummings 
model. We also discuss the influence of decoherence. The presented scheme is able to drive an initially 
separable state of two qubits into an highly entangled state suitable for quantum information processing. 

PACS. 03.67.-a Quantum Information - 03.67.Bg Entanglement production and manipulation - 42.50.Pq 
Cavity quantum electrodynamics 

1 Introduction of information. Some schemes to entangle spatially sepa- 
rated, not directly interacting, pairs of qubits via single 

Entanglement being a quantum correlation between vari- photon interference effects has been proposed [21513] . In 

ous parts of a system is required for quantum information this article we perform analytical calculations of entan- 

processing. The quantum logic gates with qubits interact- gi em ent of two distant qubits by swapping [2J. To this 

ing directly with short range interaction are not suitable end we con sider a model of two quantum two-level sys- 

for linking distant nodes. Quantum networks should be tems each independently coupled to a single mode boson 

linked with light [I] which is the best long-distant carrier field ( see pig^. We operate in the regime in which the 
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Jaynes-Cummings (J-C) Hamiltonian [5] giving two sep- 
arate qubit-field entangled states. Then we subject the 
boson field mode from each pair to a joint measurement 
(BSM). This procedure projects two formerly indepen- 
dent qubits onto an entangled state that exhibits non-local Fig- x sketch of the considere d system. Each qubit Q is placed 
quantum correlations. in its own cavity R The BSM performed on photons leaving 

In one of our recent papers we have discussed the entangle- the cavities entangles the qubits. 




2 Conditional entanglement of qubits 



ment of flux qubits using the above procedure [6] . As fol- 
lows from detailed calculations, flux qubits are usually so 

strongly coupled to the electromagnetic field modes that Let us consider two separate qubit-field subsystems (Qi?)i 



the J-C model is not adequate and one has to perform and (Q R h described by the Jaynes-Cummings Hamilto- 
numerical calculations or use the higher order approxima- man 

tions |7|. However there exists a range of qubits (atoms, u , % ( t i M 

ions [5], solid state charge qubits [5]) for which the inher- 



ent qubit-field coupling is weaker. We do not focus here 
on any of the specific examples but perform some general 
model calculations valid for systems which can be well 
described by the J-C Hamiltonian. This model is exactly 
solvable and in this paper we take advantage of it and per- 
form analytical calculations of some of the entanglement 
monotones. We derive the formulas for the linear entropy 
and show that it can be related in a simple way to the 
probabilities that can be easily measured [TO] . In the first 
part of the paper all effects of dissipation and decoherence 
are assumed to be negligible over the studied time scales 
(strong coupling limit: g > k, 7; 7, k are the decay rates of 
the qubit, field respectively). We then calculate the den- 
sity matrix for the coherently coupled qubits and take into 



Hgi (acr + + aV_) (1) 

with the coupling constant <7,-, field frequency wr 4 and 
qubit frequency ojQ i . Index i = 1 , 2 numbers the subsys- 
tems. We assume gi/uiQ i < 0.1. We have checked that in 
this regime the calculations with {1} are in agreement with 
exact numerical calculations |6l7j . The eigenstates of the 
uncoupled (gi = 0) Hamiltonian ([I]) are tensor products 
of the qubit and the field states | | n) = | |) ® \ n ) an d 
I i n ) — I I) ® \ n )\ they describe the qubit in excited | [) 
and ground | ]) states with a defined photon number n. 
The interaction term couples these states separating from 
the Hilbert space two-level subspaces S n {\ I n), | t n+1)}. 
Diagonalizing ([T]) at resonance ( lu^ = ujQ i = lur) we ob- 
tain two eigenstates 



account the influence of the decohering environment. n ^ — \ (| |, n > +| |, n + 1)) , (2) 

V2 



|,n>-|T,n + l» 
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and corresponding energies the Bell state and taking trace over the photonic degrees 

1 of freedom: 

-E± = (n + l)u R ±g. (3) 

In general n can be an arbitrary integer number but the ^ = ^ (1^^))) (9) 

Hamiltonian Q} couples only \n) and \n + 1) states. For Xo cons truct the projector we have chosen 
simplicity we assume n = that reduces S n to 



|r>B = l/V2(|01>-|10» 

{| t 1), | I 0)}. The entire system at t = is described by 

the vector state, because of the easiness of its experimental verifica- 

|-0(O)) = |i/>(0))i ® |^>(G)) 2 , (4) tion. Then the resulting qubit-qubit state reads 

where \ip(0))i describes the relevant {QR) l subsystem. |qq^ _ e -*("%+"H 3 )* r cos ( 5l £) cos(g 2 t)\ |T> 

We discuss the entanglement for two different initial states: . . . . . 

8 + S m{ gi t)sm(g 2 t)\ U)]. (10) 

hKo)> = U o)i ® 1 1 1)2 (5) 

After normalization the qubit-qubit density matrix pqq 
and is given by: 

hTO-UHh-Uflh. (6) p aa - (ID 

As the calculation procedure goes is the same way for both To quantify the strength of quantum qubit-field and qubit- 
initial states we present below a detailed analysis only for qubit correlations we calculate the linear entropy Sl [II] 
the first case ([5]) giving merely the resulting formulas for 
the second one. 

The unitary evolution of the (QR) t s generated by © leads where Pred = Ttb [pab] is the reduced density matrix, 
(for g^ ^ kir/2, k integer) to entanglement of the qubit S L = for disentangled states and reaches 0.5 for maxi- 
and field states mall y entangled states. 

Using this formula we obtain the qubit-qubit linear en- 

m)h = e~ luJRlt (cos( ff it)| | 0) -ism( 9l t)\ f 1)) , 

tropy S L 

\i;(t)) 2 = e-^ t (-ism(g 2 t)\i0)+cos(g 2 t)\n))-(7) 

S L = i- cos4 (giQ cos 4 (g 2 t) + sm 4 { gi t) sm 4 (g 2 t) 
During the evolution, the two (QR)i systems do not in- {Tr[pQq\) 

teract with each other and their state remains separable where 

|V(*)> = M*)>i |V(*))a- (8) Tr [p QQ ] = cos 2 (. 9l i) cos 2 {g 2 t) + sm 2 ( 9l t) sm 2 (g 2 t). (14) 

To entangle the qubits one needs to perform the BSM on We are also in position to calculate the qubit-qubit den- 
the field modes R\ and R 2 by projecting \tp(t)) onto one of sity matrix pqq which can be reconstructed in quantum 



(Predf (12) 
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tomography experiments |10j . and 

(IT IpqqI IT) = p. (15) (Tl \ PQQ \ IT) = (IT \pqq\ Tl) = - M2git) sin(2 ^ } 



1Tr[p QQ ] 

cos 2 (git) cos 2 (g 2 t) ^ 25 ^ 

cos 2 (g!t) cos 2 (g 2 t) + sin 2 ( gi t) sin 2 (g 2 t) Again P 1 = 0, P 4 = and the relation (J20J) is also true. 



(U\PQQ\U)=P3 (17) 

P 3 = Tt . g^M^!L : (18) 3 Results 

cos-°(5it) cos^(52*) + sm (gii)sin (g 2 i) 

/tl I I I 1 1 \ — s ' m ( 2 9it) sm(2g 2 t) 

\ I i \PQQ\ 1 1 / — U I iPQQl 1 1/ — 2Tr[pQQ\ ^ n ^ e ^ 1S ^ P ar ^ °^ section we present results obtained 

mi • . • i rrii • , r from above formulas for coherent evolution of the QR sub- 

Ine remaining matrix elements are zero. Ine signature ol ^ 

, , ,, 1,1 n systems. The influence of dissipation is discussed in the 

entanglement are tne non-diagonal matrix elements, Be- ■' 

* iv. u u-i-j.- r> ■ 1 . a jii, r 4- second part. The presented results are for the resonant 

tween tne probabilities rj, % = 1 — 4 and tne linear entropy 

there exists a simple relation case UJ< ^ i ~ URi ~ UJr ' 

Let us first assume g\ = g 2 = g; for concreteness we 

Sl = 2P 2 P 3 (20) 

take e.g. g/ujR = 0.01. Fig. [5] depicts the results for the 

Because we are working with the J-C Hamiltonian and initial state S iven by ©• At the top part we show the 

due to the projection onto \j,-) R state only two time evolution of the qubit-field linear entropy for (QR)i 

(| |T)and | tl)) from the four qubit-qubit states have finite ( Sl (qr)i ' dotted line ) and ^ R ) 2 ( S hQR) 2 ' °P en s q uares ) 

probabilities subsystems. The third curve (Sl, solid line) shows the 

For the case the system starts from the initial state qubit-qubit linear entropy as a function of time (hence- 

1^(0)) = | | 0)i ® | | 0) 2 (or |^(0)) = I T l)i ® I T 1)2) forth callcd thc ' BSM timC,) at which the BSM haS been 

and using the same procedure we obtain the following for- Performed. It shows the strength of the QQ correlations 

mula for thc linear entropy after the BSM Performed at the moment u R t. We see that 

the QR entropies overlap and their shape exhibits oscilla- 



Sl = 1 _ sin 4 (gn t ) cos 4 (g 2 t) + cos 4 (git) sin 4 (g 2 t) ^ 

(Tr [pqq]) 2 ' tions. The value of Sl is also oscillatory and depends on 



rr r 1 • it ,\ it .\ , 21 ,\ ■ 2/ ,\ the values of 5i,„„, i.e. on the degree of entanglement of 

Tr [pqq] =sm (git) cos' (g 2 t)+ cos z (git) sm (g 2 t) (22) * 6 6 

both subsystems. The maximal qubit-qubit correlations 

The corresponding probabilities are given by 

are obtained by performing the measurement when the 



P 2 = — g — g (23) QRs are maximally entangled. The Sl can be treated as 



sin 2 (ffif) cos 2 (g 2 t) 
sin 2 (gii) cos 2 (g 2 t) + cos 2 (git) sin 2 (g 2 i) 

a map of QRs entanglement onto the QQ entanglement: 

p = cos 2 (git) sin 2 (g 2 t) 

3 sin 2 (git)cos 2 (g 2 t) +cos 2 (git)sin 2 (g 2 t) the more entangled are Qi?s the more entangled are qubits 
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after the BSM. 

The QQ entanglement is also reflected in the occupation 
probabilities - it is presented in the bottom part of Fig. 
[2J The maximum of Sl, obtained for gt = kit/ 4, corre- 
sponds to equal probabilities of finding the qubits in | fj,) 
and | |t) states i.e. the qubits are then in the Bell state. 
The interesting situation arises for the initial state ©■ 



0.5 
0.4 



o 

e 0.3 

a 0.2 
3 0.1 

1 

& 0.8 

I 0.6 
ra 

■8 0.4 
£ 0.2 




3 p 1 ffi 


cfi7 P 
i/ I 1 ? 




/I / 1 


6 / \ ^ S"l QR)z 


□ 


P / \ 




4 M > 


P / \ 
fly \ 




• / \ • 

• / \ • 


• / & 


• 




7 p 3 - 

7 P 1 = P 4 


= 


/• »\ 


/• 




/ • • \ 


/ • 





100 



200 



300 



0) R t 



Fig. 2. (color online) Top: The linear entropies of qubit-field 
Sx (QH)i (dotted line), Sl (qr)i (open squares) and qubit-qubit 
Sl (solid line). Bottom: The occupation probabilities of the 
qubit-qubit state after the BSM; Qi/uor = 0.01, ujq^/ujr = 1. 

As follows from Eqs (|21j|24p for 171 = 52 = g the formulas 
for Sl , Pi and P3 reduce to 



s _ 1 / sin(2. 9 V 1 



2 \sin(2gt) J 2 



P2 = Pz 



1 / sin 



(26) 



2 Vsm(2^) / 2' 

i.e. after the BSM we always (except for gt = kir/2) get the 
qubits in the maximally entangled state. This is visible in 
Fig. |U for .92 = 0.01. For gt = kn/2 the vectors \&(t))i and 
YP{t))2 represent separable states and the state vector of 



the whole system \tp(t)) has no non-zero components along 
the direction of the Bell projector and thus the BSM is 
unsuccessful. 

Next we investigate the behavior of the system for dif- 
ferent values of gi in both (QR)i subsystems. In Figs [3] 
and H] we present the evolution of the linear entropy Sl 
as a function of the BSM time for fixed gt/uiR = 0.01 and 
varying g 2 for the initial states | J. f 1) and | J, J. 0) 
respectively. As one can see, these figures exhibit similar 
oscillations of entanglement in the almost whole range of 
parameters except for g\ = gi . In this case in Fig. [4] the 
BSM results in the maximally entangled Bell state for ev- 
ery try of measurement for which there exists the non-zero 
component of in the direction of the Bell projector. 
When gi is sufficiently far from gi the resulting entangle- 
ment is more or less regular function of the BSM time. 

Now we discuss the density matrix pqq. Its elements 
are shown in Fig[S]and Fig. [7]for the BSM time t = 0.2/is 
and for the initial states | J, f 1) and | J, J, 0) re- 
spectively. These results convincingly show signatures of 
an entangled QQ state namely the diagonal and non-zero 
off-diagonal matrix elements. We have chosen a relatively 
long BSM time so that the decoherence effects, which 
will be calculated below, are already visible. Till now 
we have presented calculations of the coherent evolution 
of the investigated system. Coupling to additional uncon- 
trollable degrees of freedom leads to various decoherence 
processes in the qubit-field evolution which results in gen- 
eral in mixed states that can be described by the master 
equation in the Markov approximation. Following [5] we 



M. Kurpas, E. Zipper: Entanglement swapping between electromagnetic field modes and matter qubits 



0.03 
0.025 
0.02 
ra 0.015 
0.01 
0,005 





100 200 300 400 500 
BSM time 



600 



Fig. 3. (color online) The dependence of the qubit-qubit linear 
entropy Sl on the BSM time for different values of the coupling 
constant Q2 for the initial state | | | 1), g\/ujR = 0.01. 
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Fig. 4. (color online) The dependence of the qubit-qubit linear 
entropy Sl on the BSM time for different values of the coupling 
constant Q2 for the initial state | J. J, 0), Qi/ur = 0.01. 

assume that the effect of environment can be included in 
terms of two independent Lindblad terms: 



/*>/,•(/)= [L H -\L n -^L K )i„ )lt U) (27) 



where the 'conservative part' is given by 



whereas the 'Lindblad dissipators' 

L m (.) = A\ n A m {-) + {-)AlA m 
— < 2A m (-)A m 



(29) 



MO 



-i\H. 



QRi 



(28) 



are expressed in terms of creation and annihilation oper- 
ators 'weighted' by suitable decoherence rates A 7 = 
and A K — <j_^/~k, m = K, 7. For concreteness we assume 
ljr = 10GHz, 7 = 1MHz, k = ZMEz. 

We have shown [5] that the BSM applied to the den- 
sity operator of mixed states is a well defined operation 
of projection and reduction which is completely positive 
and thus applicable to arbitrary density operator. Thus 
the whole procedure of swapping can be also performed if 
the decoherence effects are taken into account. 
The results of the calculations of the matrix elements pqq 
which follow from the solution of the master equation 
are presented in Fig|S] and FigJS] We see that dissipa- 
tion causes significant decrease of the non-diagonal ele- 
ments and the emergence of the additional matrix element 
(TT \pqq\ TT) — Pi 7^ that means that the two qubits 
are already not in a pure state, but are still entangled. 

Dark counts in the detectors and imperfect mode match- 
ing of photons on the beam splitter can also reduce the 
fidelity of this scheme. However it has been shown [3] that 
this kind of entangling operations is rather robust against 
such errors. Another important aspect of entanglement 
is the evolution of the state of correlated qubits. In the 
case discussed in this paper the QQ state created by the 
BSM at certain moment t does not evolve in time when 
neglecting decoherence. With decoherence taken into ac- 
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the entanglement process. However, this problem has been 
studied in some papers (see e.g. [T2]). 



Fig. 5. (color online) The qubit-qubit matrix elements at the 
BSM time t — 0.2/is for coherent evolution of QRs. The initial 
state \eOgl), UQ i = u Ri = lur, gi/ujR = 0.01, lc r = 10GHz. 




itt> \m \w 

Fig. 6. (color online) The qubit-qubit matrix elements at the 
BSM time t — 0.2/is for dissipative evolution of QRs. The 
initial state |e0<?l), coQ i = WR i = lor, gt/tOR = 0.01, lur = 
WGHz. 

count the amplitudes of the diagonal and off-diagonal ele- 
ments of pqq decrease at the expense of the appearance of 
the qubits in the ground state. In this paper we have not 
considered this problem in detail, concentrating mainly on 




Fig. 7. (color online) The qubit-qubit matrix elements at the 
BSM time t — 0.2/is for coherent evolution of QRs. The initial 
state |e0e0), LJQ i = uj Ri = ujr, gt/cjR = 0.01. 



4 Entanglement of spins encoded in quantum 
dots 

Similar considerations can also be used to entangle qubits 
encoded in the electron spin of individual quantum dots 
as recently proposed in [13j (the spin states are very long 
lived with relaxation times of order of milliseconds). In 
this approach (Fig. the establishment of spin-photon 
entanglement occurs through conditional (on the spin ori- 
entation) Faraday rotation of the incoming photon in a 
micro-cavity. Using the notation from [13) we obtain 



mn = i w t )>i ® \#ph Q (T)h 



(30) 
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Fig. 8. (color online) The qubit-qubit matrix elements at the 
BSM time t — 0.2/is for dissipative evolution of QRs. The 
initial state |eOeO), UQ i = u)^ = cjr, gi/tUR = 0.01. 
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Fig. 9. Entanglement swapping procedure for photons and 
spins in quantum dots. The spin-photon interaction produces 
a conditional single-photon Faraday rotation. For the rotation 
angle n/4 the spin-photon state is maximally entangled. The 
BSM performed on such two photons produces entangled spins. 

|<WT)>i = VV2(l\>ilT>i + l/">iU>i) 

\%h Q {T)) 2 = l/x/2 (| \) 2 | t> 2 + I /> 2 | i) 2 ) , (31) 

where T is the interaction time between the electron spin 
and the photon in the micro-cavity, | \,) and | y) the pho- 
ton states with a linear polarization rotated by — 7r/4, +7r/4 



respectively with respect to the state | <->) of linear po- 
larization in the x direction. The BSM on the photons 
outgoing from the two micro-cavities conditionally leads 
to entangled qubit states 

\#QQ) = Tr vh l\*-) vhph {*-\*){9\) 

= -l/V2(|t>iU>2-U>i|t>2) (32) 

where 

\#-) ph = l/x/20 \),\ Sh - | \) 2 ) (33) 

In this case the polarization degrees of freedom make up 
the photonic qubit. This scheme provides a link between 
spintronic and photonic systems and can combine the ad- 
vantages of both spintronic and photonic quantum infor- 
mation processing. 

5 Conclusions 

Cavity quantum electrodynamics with individually address- 
able qubits (atoms, trapped ions, charge, flux or spin solid 
state qubits) is expected to provide a toolbox for quantum 
computing. The strong qubit-field coupling achievable in 
a high-finesse cavity can be accurately described by the 
Jaynes-Cummings model if g/ujQ < 0.1. This interaction is 
coherent permitting the transfer of quantum information 
between the qubit and the electromagnetic field modes. 
Photons are natural candidates as carriers of quantum in- 
formation because they are highly coherent and can me- 
diate interaction between different objects. 
We have considered two remotely located stationary qubits 
each of which becomes entangled with the respective elec- 
tromagnetic field modes. Such qubit-field entanglement 
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has been demonstrated e.g. for charge qubits [9] and for References 



atoms j8] . When the electromagnetic field modes are com- 
bined on a beam splitter their appropriate coincidence 
measurements ensure the entanglement of the two qubits. 
This procedure can be used both for pure and mixed qubit- 
field states and the results depend on the initial state of 
the system. The created entanglement survives in the pres- 
ence of dissipation processes with the realistic decoherence 
rates of both qubits and fields. 

The fact that the generation of the entangled state of 
qubits is conditional on the BSM on photons to some ex- 
tent limits the application of this technique. However this 
is not an obstacle to what seems to be the most promis- 
ing application - the entanglement of two distant qubits 
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